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ABSTRACT 

Formation of primordial black holes (PBHs) on astrophysical mass scales is a natural 
consequence of inflationary cosmology if the primordial perturbation spectrum has a 
large and negative running of the spectral index as observationally suggested today, 
because double inflation is required to explain it and fluctuations on some astrophysical 
scales are enhanced in the field oscillation regime in between. It is argued that PBHs 
thus produced can serve as intermediate-mass black holes (IMBHs) which act as the 
observed ultraluminous X-ray sources (ULXs) by choosing appropriate values of the 
model parameters in their natural ranges. Our scenario can be observationally tested 
in near future because the mass of PBHs is uniquely determined once we specify the 
values of the amplitude of the curvature perturbation, spectral index and its running 
on large scales. 
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1 INTRODUCTION 

Ultraluminous X-ray sources (ULXs; Makishima et al. 
2000) are characterized by their luminosity greater than 
~ 10''® erg s~^ (well above the Eddington luminosity of a 
neutron star) and off-nuclear location in nearby galaxies 
(Fabbiano 1989; Colbert & Mushotzky 1999). In addition, 
X-ray variability on various timescales are found for some 
ULXs (Matsumoto & Tsuru 1999; Ptak & Griffiths 1999; 
Kaaret et al. 2001; Kubota et al. 2001). Luminosity, time 
variability and X-ray spectra of ULXs indicate that ULXs 
are acc reting black holes, r ather than young supernova rem- 
nants l|Kaaret et al.l [2OOII ). Although association of ULXs 
with active star-forming regions are clearly shown (Mat- 
sushita et al. 2000; Zezas et al. 2002), the physical reason 
behind the association, and moreover the origin of ULXs, 
are still unclear. 

The number density of ULXs can be roughly esti- 
mated as follows. Late-type/starburst galaxies hosting nu- 
merous compact X-ray sources (mostly high-mass X-ray 
binaries) within about 30Mpc from the M ilky Way are 
listed in (|Grimm. Gilfanov fc Sunvaevj|2003l ). The number 
of these X-ray sources with each luminosity greater than 



2x lO'^^^erg s~^ is about 100. Grimm et al. also show that the 
luminosity function N{L) of compact X-ray sources, which 
is defined as a number of X-ray sources with luminosity 
greater than L, is universal for galaxies if it is normalized 
proportionally to the star-forming rate (SFR). The normal- 
ized luminosity function tells that among L > 10^^'^erg s~^ 
sources about one-third have L > 10^^ erg s~^. Hence, the 
number of ULXs with L > 10^® erg s"^ wiU be ~30 within 
30Mpc. Thus, the number density of ULXs is estimated to 
be 10"^'^Mpc-^ 

In spite of observational evidences, theoretical expla- 
nations of ULXs are still in dispute. The most challenging 
fact is the high luminosity greater than the Eddington lu- 
minosity of a stellar-mass black hole with mass ~ IOMq. 
For example, it is argued that sources of ULXs can be stan- 
dard stellar-mass black holes with jets or relativistic beam- 
ing l|Koerding. Falcke fc Markofj [2OO2I ). Another possibib 
ity is stellar-mass black holes radiating at super-Eddington 
luminosities due to efficient photon leakage from accre- 
tion disks (e.g., Begelman 2002; Meyer 2004) and/or due 
to super-Eddington accretion rates ("slim accretion disks"; 
Abramowicz et al. 1988). Detailed computations of slim ac- 
cretion disks (Kawaguchi 2003) indeed account for observed 
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X-ray spectra of ULXs (Okajima, Ebisawa & Kawaguchi 
2006; Vierdayanti et al. 2006; Ebisawa & Kawaguchi 2006; 
Foschini et al. 2006). 

ULXs can also be explained as sub-Eddington accre- 
tors by assuming black holes heavier than stellar-mass black 
holes (see Miller & Colbert 2004 and reference therein). 
Moreover, ULXs exhibit quasi-periodic oscillations (QPOs) 
at frequencies lower than the QPO frequencies of normal 
black hole binaries by a couple of orders (Strohmayer & 
Mushotzky 2003; Strohmayer et al. 2007). A straightforward 
interpretation of the longer timescales is a heavier black 
hole mass than IOMq. Hence, intermediate-mass black holes 
(IMBIfs) with mass ~ are considered to be one of 

the best candidates of ULXs. 

The most challenging problem with such IMBHs is their 
formation mechanism. There are several possible mecha- 
nisms discussed so far. For instance, detailed evolution mod- 
els of stellar binaries show that the generation rate of IMBHs 
is very small (Madhusudhan et al. 2006). 

A promising possibility is that IMBHs could be rem- 
nants of Population III stars (e.g., Schneider et al. 2002). 
It is suggested that zero- met alicity stars form with masses 
1O^"^M0 (e.g., Omukai et al. 2005). If the initial mass of 
a star is sufficiently large, it can collapse directly to an 
IMBH. These IMBHs can explain the observ ed nurnber o f 
ULXs if they are generated in galactic disks (Krolikl l20o3 ). 
It is also discussed that IMBHs produced in galactic halos 
can account for obse rved ULXs if Popula tion III star forma- 
tion is very efficient (iMii fc Totanill2005^ . This possibility is 
suggested by the e xcess of the cosmic near-infrared back- 
ground radiation (I Wright fc Rees3 l2000l : ICambresv et al.l 
I2OOII : iMatsumoto et al.l l2005h . Note, however, that with 
the typical baryonic fraction of Pop-Ill stars, / ~ 10~^ 
jMadau fc Reesl l2001^ ■ the abundance of remnants BHs, 
Qbh = f^b ~ 4 X 10~^, would be too small to account 
for IMBHs as we will see below (§4). On the other hand, if 
we assume a larger baryonic fraction, / > 1 0~^, then the 
Univ erse may have been reionised too early (jPaigne et al.l 
I2OO6I ). 

Black holes can be produced via the collapse of over- 
dense region of t he ea,rly un iverse with large initial curva- 
ture fluctuation (|Carilll975l ). These black holes are called 
primordial black holes (PBHs). If the power spectrum of 
initial curvature perturbation has strong peak at a specific 
scale, collapse takes place at the epoch when this scale en- 
ters the horizon. Typical mass of resultant PBHs can be 
determined by the horizon mass in this epoch. Formation of 
PBHs in inflatio n ary cosmology is discussed in (|lvanovll99 7 ; 
Yokovama"l997; Kawasaki ct al. 1998; Y okoyamalll998al lb: 
Kawasaki fc YanagidaT999, : , Yamaguchi . 200 ll ) . 



One exotic possibility is the formation of IMBHs dur- 
ing^ Quantum Chromo Dynami cs (QCD) pha se transition 
(fjedamzik 1997, 1998; Jedamz ik fc Niemeveij [19991 . if the 
transition is flrst order. Since the equation of state of the 
universe is relatively soft in this epoch, collapse of overdense 
region takes place easily. However, the horizon mass at this 
epoch is only 0{1)Mq rather than lO^'^M©. 

In this paper, we present a model of inflation 

I Yamaguchi fc Yokovamal2004l : [Kawasaki et al ]|2006l ) which 
can produce significant amount of PBHs in astrophysically 
interesting mass scales, to show that the origin of IMBHs 
could be such PBHs. In this scenario, the spectrum of ini- 



tial curvature pert urbation has strong peak due to the para- 
metr ic resonance l|Kofman et al ] |1994| . Il997l : IShtanw et al.l 
After infiation, the oscillation of the infiaton con- 
densate can result in oscillating effective mass of another 
scalar field and/or the infiaton itself. This oscillating ef- 
fective mass excites large fiuctuations of that scalar field 
for specific modes. If this scalar field contributes to the 
energy density of the universe, excited fiuctuation results 
in large curvature perturbation for this specific modes. If 
this curvature perturbation is suffici ently large, PBHs are 
produced via g ravitational collapse (| Green fc MalikI [2OO1I : 
iBassett fc Tsui ikawa 2001). In general, characteristic scales 
of fiuctuations excited by parametric resonance corresponds 
to very short wavelength compared to present observable 
scale. However, because another infiation is required after 
the parametric resonance in our model, these large fiuctu- 
ations are naturally expanded to cosmologically and/or as- 
tronomically relevant scales. 

Interestingly, this scenario is related to the running 
spectral index of the universe. Although it is not compelling, 
the result of the Wilkinson Microwave Arusotropy Probe 
(WMAP) observations (Komatsu ct.al. 2008) indicates that 
the running of spectral index may be large and negative. 
These large running of the spectral index cannot be ex- 
plained unless two or rnore successive inflation s are assumed 
iKawasaki et al.ll2003l : [Easther fc Peiri^ '2006V Therefore, if 
a large running spectral index is confirmed by forthcom- 
ing observations of cosmic microwave background (CMB), 
strongly peaked initial curvature perturbation from multi- 
ple infiation can be an interesting candidate of the formation 
mechanism of IMBH. 

This paper is organized as follows. In Section 2, we ex- 
plain the scenario of infiation and how the strongly peaked 
initial fiuctuation is produced. The calculation of the mass 
distribution of PBHs is given in Section 3. In Section 4, we 
estimate the amount of IMBHs produced via the collapse 
of over dense region. In Section 5, we discuss the relation 
between the running spectral index and the typical mass of 
PBHs. The possibility of explaining other astronomically in- 
teresting PBHs is also discussed. Finally, we summarise the 
results. 



2 GENERATION OF PEAKED POWER 
SPECTRUM OF INITIAL CURVATURE 
PERTURBATION VIA A DOUBLE 
INFLATION MODEL 

2.1 Smooth hybrid new inflation model 

The power spectrum of initial curvature perturbation is 
needed to be strongly peaked in order to result in the for- 
mation of relevant amount of PBHs via collapse of over- 
dense region. Initial curvature perturbation is determined 
by the detail of inflation. Here we consider the smooth 
hybrid new inflation model in supergravity proposed by 
lYamaguchi fc Yo kovama! (|2004l ). It was first considered to 
account for the running spectral index suggested by WMAP. 
The spectral index n of the initial curvature perturbation 
Vtz is suggested to be dependent on the comoving wavenum- 
ber k, according to the result of WMAP observation. The 
best-fit values of the amplitude Viz, spectral index n and its 
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running a = dn/d\nk at the pivot scale fco ~ 0.002Mpc ^ 
are 



Vn = (2.40 ± 0.11) X 10'^ n = l.OSlto.oss , 
a = -0.037 ±0.028 



(1) 



by the WMAP 5 year result only (|Komatsu et.al.l |2008| ). 
This suggests the existence of large running from n > 1 at 
long wavelength scale to n < 1 at short wavelength scale 
within the range WMAP observation can probe, namely, 
IMpc < fc-^ < lO^Mpc. 

In our model, two inflatons are introduced and 
two stages of i nflatio n takes place in succession. (See 
iKawasaki et all (j2006l ) for the more detailed discussions 
based on supergravity.) At first, the smooth hybrid infla- 
tion occurs. This inflation is characterized by the following 
effective potential of an infiaton real scalar field a, 



V{a) 



1 - 



8^i' 2 



2 ifiMf 
27 cr4 



+ 



8M^ 



for cr > {fiM) : 



for a < (^iM) 2 



Here parameters with dimension of mass /i and M indicate 
the energy scale of smooth hybrid infiation and the cut- 
off scale of the underlying theory, respectively. Here Ma is 
the reduced Planck mass, Mg = 2.4 X lO^^GeV. Typically, 
parameters /i and M satisfy i-i <^ M < Mg- The infiaton cr 
slowly rolls from cr ^ (/lA/)^''^ to cr = 0. Here you should 
notice that the observational quantities Viz, n, and a at 
the pivot scale fco of primordial fiuctuations are completely 
determined by the model parameters n, M, and cr = ctq at 
the moment the scale fco exits the horizon. The duration 
of the smooth hybrid infiation after the scale fco exits the 
horizon is also determined once crini, /i, and M are specified. 

Smooth hybrid infiation model is a variant of hybrid in- 
fiation model. In both models infiation is driven by a false 
vacuum energy of some symmetry-breaking field. In hybrid 
infiation, the restored symmetry is broken only at the end of 
infiation resulting in formation of topological defects which 
could be cosmologically harmful. On the other hand, dur- 
ing smooth hybrid infiation, by virtue of non-renormalizable 
terms, symmetry remains broken and there is no defects 
formation at the end of infiation. These non-renormalizable 
terms can also serve to produce appreciable negative running 
spectral index. Furthermore, while it has been shown that 
hybri d infiation requires severe fine-tuning of initial condi- 
tion, iTetradis N.lll998l : iMendes L. E]|2000h smooth hybrid 
infiation is free from such a fine tuning and can be realized 
naturally. 

The most important feature of this infiation model is 
that the largelly runiung spectral index can be realized. This 
potential has positive curvature for larger cr and negative 
curvature for smaller cr. Since the spectral index of the ini- 
tial curvature perturbation depends on the gradient of the 
potential, this infiation results in scale-dependent spectral 
index. Large wavelength modes of the initial curvature per- 
turbation are produced while the infiaton takes a large value. 
Therefore, the spectral index for these modes is larger than 
unity. On the other hand, short wavelength modes are pro- 
duced while the infiaton takes a smaller value. Thus, the 
spectral index for these modes is smaller than unity. Con- 
sequently, negative running of the spectral index suggested 



by the WMAP 5 year result can be realized. For example, a 
parameter choice (equations are corrected) 



= 2.1 X 10"'^ Mg, M = 1.3Mg, ao = 0.227Mg 



gives 

Viz = 2.40 X 10" 
a = -0.033, 



n = 1.040, 



(3) 



(4) 



which is in the range ([T]). However, since large running needs 
large variation of the gradient of the potential, the evolution 
of infiaton must be very fast in order to account for the large 
running of the spectral index, which implies that smooth 
hybrid infiation does not last so long. Thus, another infiation 
is necessary to push the relevant scales to cosmologically 
observable scale. 

In this model, new infiation is considered as another 
infiation to compensate for the duration of the infiationary 
epoch. Moreover, it naturally leads to low reheating tem- 
(2) perature to avoid the overproduction of the gravitinos. The 
effective potential of the second infiaton field (j) during this 
new infiation is given by 



2 Ml 



2 Ml 16 MA ' 



(5) 



where v determines the scale of new infiation, which is as- 
sumed to be n <C Parameters c < 1 and g < 1 deter- 
mine the shape of the potential V{(j)). In addition, the sec- 
ond infiaton field <j) is coupled to a through gravitationally 
suppressed interactions in supergravity even if no direct cou- 
pling is introduced. Since the contribution from this coupling 
dominates the potential of 0, 4> is trapped to a value close 
to, but different from <j!> = 0, while cr has a non-vanishing 
value. After the smooth hybrid inflation, a oscillates around 
cr = with its amplitude decreasing. Eventually, the con- 
tribution from V{(l)) becomes dominant, and then begins 
to roll slowly from the vicinity of </> = to the minimum 
(f) ~ {Av^ Mg/ gY^'^ ■ The duration of new inflation is de- 
termined by parameters v, g, and c. Note that fluctuation 
generated during the new inflation is not constrained by ob- 
servations unless it is extremely large, because length scale 
of such fluctuation is too small to be cosmologically observ- 
able. The exception is th e strong fluctuation g enerated via 
the parametric resonance (|Kawasaki et al.ll2006l ). which may 
serve as a seed of PBHs, as discussed below. 



2.2 Generation of strong fluctuation via the 
parametric resonance 

Because of the self-coupling of a, oscillating condensate of 
a gives oscillating contribution to its own effective mass. 
Since this self-coupling is large, the amplitude of this os- 
cillating contribution is larger than its mass at the origin 
rria ~ {Sfi'^/MY^^. Hence some speciflc modes of fluctuation 
of a are strongly amplifled via the parametric resonance. 
In our case, it is determined by the evolution equation of 
Fourier modes of fluctuation at , which can be approximated 



CTfe + [A- 2gcos(2z)] at = 0. 



(6) 



Here A is a function of k/{ama), where a is the scale fac- 
tor, g is a time-dependent parameter proportional to the 
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self-coupling and the amplitude of the oscillation of the 
zero-mode of a. The prime represents derivative with re- 
spect to the time variable z defined by the proper time 
as 2z — mat — n/2. Equation ^ has the shape of Math- 
ieu equation. The solution of this equation have infinite 
sequence of instability bands determined by the value of 
A. In these instability bands, the solution grows exponen- 
tially for g > 1. Hence, the modes k satisfying the condi- 
tion k/a — ma X (const.) grows exponentially. According to 
numerical calculation, the strongest instability appears at 
around k/a = x 0(10"^). Note that the comoving scale 
k in this instability band is inside the horizon at the begin- 
ning of the new inflation. Eventually q decreases to g ^ 1 
and the parametric resonance ceases. Later, fluctuations (j>k 
with the same range of k get large amplitudes via linear 
coupling with the amplifled fluctuation at- These fluctua- 
tions (f>k determine the curvature perturbation during the 
new inflation. Thus, the curvature perturbation has a large 
amplitude at a speciflc range of k. 

We can see that the position of the strong peak kp is 
approximately determined by the factor {/j-M) and the du- 
ration of the smooth hybrid inflation. The peak position kp 
satisfies 



: Cma 



(7) 



where C ~ 0(10"^) is determined by the detail of the para- 
metric resonance and is the scale factor at the end of the 
smooth hybrid infiation. On the other hand, the scale fac- 
tor ao at the moment the pivot scale exits the horizon and 
the Hubble parameter during the smooth hybrid infiation. 



Hs ~ fM^ /(VSMg), satisfies 
ko 



aoHs 



Therefore, we get 
ko ao 



Mg 



+ lnC+1.6. 



(8) 



(9) 



The shorter smooth hybrid inflation we assume, the larger- 
mass PBHs are produced. 

Note that the typical mass of PBHs is determined once 
we specifled Viz, n and a. In our model, only the smooth- 
hybrid inflation is responsible to the curvature perturbation 
of cosmologically interesting scales. Since free parameters of 
our smooth-hybrid inflation model are p, M and the value of 
inflaton ctq at the moment the pivot scale ko crosses outside 
the horizon, cosmological parameters Viz, n and a depend 
only on parameters /i, M and ao. Therefore, in order to 
reproduce specific values of V-n, n and a, we look for appro- 
priate values of /i, M and ao. Once we find them, a^jao is 
determined by the dynamics of the inflaton, and then kp/ke 
is specified by the Eq.([9ll. 

The amount of PBHs depends on the height of the peak, 
which is determined by the balance of the efficiencies of the 
parametric resonance and the subsequent decay of a into 
lighter particles. If this decay rate is sufficiently large, it 
makes the parametric resonance weaker. It is extremely dif- 
ficult to estimate the magnitude of this strong peak analyt- 
ically because of the non-perturbative nature of the para- 
metric resonance. 
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Figure 1. The spectrum leading to the formation of IMBHs from 
the smooth hybrid new infiation model. 



2.3 Spectrum of the curvature perturbation for 
the formation of IMBHs 

In Fig[TJ we show an example of the spectrum of the initial 
curvature perturbation which results in the formation of sig- 
nificant amount of IMBHs. In order to give a example of gen- 
eration of IO^A/q PBHs, we choose ^ = 1.6 x 10"^Mg, M = 
O.SIMg, V = /i/4, c — 0.1, and g — 10~^ for parameters of 
the smooth hybrid new infiation as an example. The position 
of the peak is at kp — 4.4 x lO^Mpc"^, thus the typical mass 
of resultant black holes is estimated as Mpbh = 108Mq by 
the horizon mass at the moment the scale kp enters the hori- 
zon. Unfortunately, it is difficult to estimate the range of 
these parameters which results Mpbh ~ IOOMq. As a hint 
for this estimation, we tried various parameters for and 
M. For example, /x = 1.7 x 10"^Mg and M = 0.86Mg 
gives MpBH IOOOM0, while n = 1.5 x 10"^A'/g and 
M = 0.72A'/g gives AfppH =i WMq, both reproducing the 
cosmological parameters consistent with ((l|. We also as- 
sumed the decay rate of cr as P = 7.65 x lO^^Afc. Modes 
of curvature perturbation with k < lO^Mpc"^ are gen- 
erated during the smooth hybrid infiation. This spectrum 
gives Vtz = 2.4 x 10"^, n = 1.067, and a = -0.014 at 
ko = 0.002Mpc~^. A large and negative running spectral in- 
dex is a characteristic feature of formation of IMBHs via this 
scenario, which can be tested by forthcoming observations 
of CMB. 



3 FORMATION OF PRIMORDIAL BLACK 
HOLES VIA COLLAPSE OF OVERDENSE 
REGION 

PBHs were considered to be produced via gravitational 
collapse if the density perturbation S = 5p/p within 
a patch of the universe is larger than some critical 
value Sc at the moment the Hu bble radiu s becomes 
smaller than the size of the patch (ICar 3I1975I). Later, a 
critical phenomenon was observed (|Niemever JedamzikI 
1998) in the formation of PBHs via gravitational col- 
lapse of radiation fiuid in the Friedma nn-Robertson- 
Walker background vans fc ColemanlUggj ). According to 
(|Niemever fc Jedam zik 1998) , the mass of PBH Mpbh pro- 
duced via collapse of a overdense region with the density 
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perturbation 5 has a scaling relation 



Mp 



(10) 



where Mh is the horizon mass at the moment this region 
enters the horizon. The index 7 and the critical value 5c 
are universal constants, which are estimated by numerical 
simulations as 7 ~ 0.35 and 5c — 0.67, where the latter 
is quote d from results of numerica l simulation of critical 
coUapse dNiemever fc Jedamziklll999h . 

Given strongly peaked spectrum of curvature pertur- 
bation, the mass spec trum of PBHs is estimated as follows 
l|Yokov ama 1998c: B arrau et al.|[2003h . Assuming that the 
initial curvature perturbation has a Gaussian probability 
distribution and that the curvature perturbation is strongly 
peaked, the differential mass spectrum of the PBHs at the 
time of formation can be given by. 



l{tc 



d In A/pBH 



1 + i 

7 



X exp 



\ Mh J 



(11) 



The parameter e is defined as the ratio of the peak of differ- 
ential mass spectrum Afmax to the horizon mass Mr- The 
function P{Mh) is the probability that the relevant mass 
scale has an above-threshold amplitude of fluctuations to 
collapse as it enters the Hubble radius, namely. 



(12) 



where p{5,tc) is the probability distribution function of the 
density perturbation at the moment the scale kp enters the 
particle horizon, which is assumed to be Gaussian. The vari- 
ance of the density perturbation auitc) smoothed over the 
comoving length scale can be estimated by 



o"H(ic 



16 
81 



k \ dk 



(13) 



Here we employed the top-hat window function 



(14) 



Note that the power spectrum of the density perturbation 
Vs is given by 

4 



Vs = 



16 



•P7^, 



(15) 



where T{k/kp,tc) is the transfer function which is given by 



k 



T ^Ac = ^ 



kp ( V^kp . / k 
— ' sm I 



fc2 \ it 



V\/3fcy 



( k 



\\l?)kp 



(16) 



Finally, the differential mass spectrum of PBHs at 
present can be given by 



dQ.Y 



\(tc 



d In MpBH d In Mpbh 



exp 



iwHdt 



(17) 



unless PBHs are dominant constituent of dark matter, where 
w is the effective equation-of-state parameter of the total 
cosmic energy density including contribution from dark en- 
ergy. In Fig[2] we show the differential mass spectrum (|17p 
calculated for the spectrum shown in Fig[T| with e = 1. 



a a 




Figure 2. The differential mass spectrum II17I I calculated for the 
spectrum shown in Fig[l] 



We took contribution from dark energy into account. We 
can see that the mass distribution of PBHs is well peaked 
around Mpbh = lOSM©. Therefore, we can conclude that 
only IMBHs are produced from this spectrum. 



4 AMOUNT OF IMBHS 

The observational constraint on the total amount of IMBHs 
can be estimated from the observed number density of ULXs 
IQ-^-^Mpc"^, as follows. Assuming that ULXs are IMBHs 
with typical mass 10^ Mq, this number density means that 
the contribution of IMBHs observed as ULXs to the density 



parameter of the universe Hulxbh is 



10" 



Agol and 



Kamionkowski (2001) estimate the probability for floating 
black holes to acquire high enough accretion rates from in- 
terstellar medium as ~ 10~^. Furthermore, if we consider 
the enhancement of dark matter (and IMBHs) in the galac- 
tic disk region properly, the typical mass fraction of dark 
matter (over the entire dark matter halo) within the galac- 
tic disk region would be ~ 10"^ (see Mil & Totani 2005). 
By combining the above three points, we can then estimate 
the total abundance of the IMBHs in the universe as 



f^ULXBH X (10~^ X 10"^)"' = 10" 



(18) 



We note that estimations adopted above would be valid only 
up to the order of magnitude. 

We here point out that JIimbh quoted satisfies the con- 
straints from CMB observations (see Ricotti, Ostriker & 
Mack 2007). Namely, if there are too many PBHs, a frac- 
tion of these black holes shining by gas accretion in the early 
universe at 2 ~ 100 — 1000 (where gas density is quite high) 
must have distorted the CMB spectrum. Then, the CMB 
spectral distortions (based on Far-Infrared Absolute Spec- 
trophotometer (FIRAS) data) provide an upper limit for 
allowed f2iMBH to be 10~*'^ for the case of the PBH mass 
AfpBH = lO^Afe (Ricotti et al. 2007). 

We calculated the abundance of the IMBHs from the 
differential mass spectrum Eq. (|17p for the same parameter 
set as we chose for Fig[l]and various decay rates F. The re- 
sult is summarised in FigO which shows that the abundance 
of IMBHs is sensitively dependent on the decay rate F. Ac- 
cordingly, the range 7.65 x IQ'^Ma < F < 7.74 x Ma 
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r [Mg] 

Figure 3. The relation between the abundance of IMBHs Cimbh 
and the decay rate F. 
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spectral index n 

Figure 4. Contour-lines of the typical mass of PBHs on the n-a 
plane. 



results in the abundance lO"'^'^ < r2iMBH < 10"^'^. Thus, 
this scenario needs tuning of the decay rate in order to ac- 
count for the abundance of ULXs. This is simply due to the 
fact that the abundance of PBHs is exponentially sensitive 
to the peak amplitude of fluctuation, and all the other mod- 
els of PBH formation requires the same level of fine tuning. 
Note that the decay rate of a is expected to be given by 
r ~ Ng^rria/ (Stt) where g^r is some gauge coupling constant 
included in near-Planck scale physics and A'' is the number 
of decay channels. Under the parameter set we have cho- 
sen, the above constraint on F is rewritten by Ng^ ~ 10~^, 
which can be satisfied by natural values of A'^ and gcr- Here 
one should notice that F has nothing to do with the reheat- 
ing temperature after the final infiation, which is determined 
by the decay rate of <f). 



5 SPECTRAL INDEX AND TYPICAL MASS 
OF PBHS 

Here let us consider the relation between the typical mass of 
PBHs and cosmological parameters, spectral index and its 
running. As we discussed in the Section 2, V-r, n and a at the 
pivot scale ko — 0.002Mpc~^ are determined once we choose 
parameters M and the value of a at the moment the scale 
ko exits the horizon. This choice simultaneously determines 
the duration of the smooth hybrid infiation ln(ae/ao). Thus, 
if we fix the amplitude Vtz to be the best-fit value of the 
result of WMAP3, we can compute the mass of PBH for each 
combination of (n,a), up to an uncertainty in the detail of 
parametric resonance. A brief explanation of the method we 
used is given in Appendix. 

In Fig|4j we show the relatiorQ between typical mass of 
PBHs and the cosmological parameters n and a. We can see 
that larger mass PBHs requires larger running of the spectral 
index, while lower mass PBHs can be produced with negli- 
gible running. This can be understood by the fact that the 

^ Cosmological parameters n and a of the spectrum shown in 
Fig[l]and the typical mass of PBHs 108Mq is deviated from this 
result. The main reason is that v is not negligibly small compared 
with II. Since long oscillatory phase requires extremely long nu- 
merical calculation, we chose large v for our calculation. 



short smooth hybrid infiation is realized for steep potential, 
which means large running. If we choose some other appro- 
priate values of model parameters, this scenario can alterna- 
tively explain one of other astronomically interesting PBHs, 
such as 0.1 — IMq BHs considered to be a candidate of mas - 
sive compact halo objects (MACHOs) (|Alcock et all |2000| ) 
or 10^^ g BHs which are p ossible sources of ul tra high energy 
cosmic rays (UHECRs) l|Takeda et al.l 120031 ') . ahhough the 
latter possibility is probably no longer int eresting given the 
recent Auger data (| Abraham et al.l l2007l ). The former re- 
quires Q ~ —0.01, while the latter requires negligibly small 
a. For the latter case, best-fit values of Viz and n is given 
by 

Ptz = (2.41 ± 0.11) X 10~^ n = O.gesto.olg (19) 

from the result of WMAP 5 year only, assuming a = 0, 
which can also be realised in the present scenario. 

6 SUMMARY AND DISCUSSION 

We considered the formation of IMBHs as a candidate of 
ULXs via gravitational collapse of overdense region of the 
curvature perturbation whose spectrum is strongly peaked. 
This initial fiuctuation is produced by a double inflation 
model, the smooth hybrid new inflation. The peak is origi- 
nated from the amplification of specific modes of fiuctuation 
via the parametric resonance. This infiation model can also 
explain the observed running spectral index. We found the 
relation between typical mass of PBHs and the cosmological 
parameters n and a. This result shows that the formation 
of IMBHs in this scenario requires significant running of the 
spectral index. The amount of IMBHs depends on the height 
of the peak, which is determined by the detail of paramet- 
ric resonance and the decay rate of the infiaton. The decay 
rate must be some specific and fine-tuned value for required 
abundance of IMBHs. The level of required fine-tuning is 
no different from other models of PBH formation. What is 
important here is that the requirement can be satisfied with 
natural values of the parameters. This scenario can alterna- 
tively explain one of other PBHs, such as 0.1 — IMq BHs 
considered to be candidates of MACHOs or lO^^g BHs which 
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are possible sources of UHECRs with some other appropri- 
ate choices of model parameters. 

The production of strongly peaked initial curvature 
perturbation via the parametric resonance can take place 
in other do uble inflation models with intermediate oscilla- 
tory phases (|Kawasaki et al.ll2003l : lYamaguchi fc Yokovamal 
l2003l i. Since these double inflation models are required to 
explain the possible large running of the spectral index, ob- 
servational confirmation of large and negative running spec- 
tral index may be regarded as a hint of IMBHs. 



APPENDIX 

Here we will explain how Fi g|4]was determined. We used an - 
alytic estimations given in I'Yamagu chi fc Yokovamall2004l ). 
because it is sufficient to see approximate dependence of 

MpBH on (n, a). 

According to l|Yamaguchi fc Yokovamal 120041 1 . the dy- 
namics of inflaton a has characteristic values of a. Namely, 
ffd, where the dominant contribution to the dynamics 
changes from the second term in Eq.(12| to the third term, 
and (Jc, where the slow-roll condition is violated and the 
smooth hybrid inflation ends. These values are estimated 
by followi ng equations. Note that we con sider only the case 
m = 2 in l|Yamaguchi fc Yokovamall2004l ). 



Wa) 
Wo) 



The number of e-foldings Nd from a — to a — ac is 



Nd 



+ 



Wo) 



(20) 
(21) 

(22) 



Now, if we choose the value of /xM and require the typ- 
ical mass of resultant PBHs Mpbh, A^e = ln(ae/ao) is deter- 
mined by Eq.(|9}. Then, the value of ctq, which is the value of 
a at the moment the pivot scale ko crosses outside the hori- 
zon, is determined by tracing back the dynamics of a from 
the end of smooth hybrid inflation a = ac. If > Nj, it 
means that ao > (Jd- Therefore, 



^0 



+ 



6 



On the other hand, if < Nf, then ao < (Jd and 



5 ,r 

6 



1 

3 



Wg) 



(23) 



(24) 



Thus, we can determine ctq if we specify (/^M) and Mpbh 
Because Viz is given by 



V37r 



3 / O-rf 

era 



I 3 



(25) 



parameters /i, and then M, can be calculated once we deter- 
mined ao- Cosmological parameters n and a are calculated 
via slow-roll parameters. 

Consequently, if we choose Mpbh, we can draw a curve 
on (/i, M) plane which results the typical mass of PBHs 



MpBH. Deriving the corresponding values of {n,a)\c=aQ, 
can draw Fig|4] 
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